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We calculate several neutron star properties, for static and/or rotating stars, using equations 
of state based on different microscopic models. These include our Dirac-Brueckner-Hartree-Fock 
model and others derived from the non-relativistic Brueckner-Hartree-Fock approach implemented 
with microscopic three-body forces. The model dependence is discussed. 
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I. INTRODUCTION 

The determination of the equation of state (EoS) of 
highly compressed and/or hot matter is one of the most 
complex problems in nuclear physics. In terrestrial labo- 
ratories, relativistic heavy-ion collisions are the best tool 
to produce hot and dense matter. Concerning astrophys- 
ical systems, neutron stars are known to be stable config- 
urations containing the most dense form of matter found 
in the universe. They are therefore unique laboratories to 
study the properties of highly compressed (cold) matter. 
Furthermore, the possibility of studying the structure of 
neutron stars via gravitational waves jlj makes these ex- 
otic objects even more exciting. 

More generally, the study of neutron-rich systems 
(from the lowest to the highest densities) has widespread 
impact, reaching from the physics of exotic nuclei to nu- 
clear astrophysics. With the Facility for Rare Isotope 
Beams (FRIB) recently approved for design and con- 
struction at Michigan State University, these studies be- 
come particularly important and timely. Partnership be- 
tween nuclear physics and astrophysics will play a crucial 
role in advancing our knowledge of neutron-rich matter 
and its equation of state. 

Clearly, a rich (on-going or future) experimental pro- 
gram must be accompanied by parallel theoretical ef- 
fort. Present nuclear matter calculations are typically 
performed either within the mean field approach or the 
microscopic one. Mean field models, relativistic and non- 
relativistic, are popular 0, Q, but in our opinion ab ini- 
tio approaches are best in order to get deeper physical 
insight. By ab initio, we mean that the starting point is 
a realistic two-body potential, possibly complemented by 
three-body forces. 

The purpose of this paper is to compare several neu- 
tron star properties predicted by different microscopic 
models. As constraints promise to become more strin- 
gent, it is important to understand and compare how 
the nature of the various predictions is related to the 
features of each model. In microscopic approaches, the 
tight connection with the underlying two-body potential 
will then facilitate the physical understanding in terms of 
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the characteristics of the nuclear force and its behavior 
in the medium. 

In the next section, we review the main aspects of the 
ab initio approach. The main differences between our ap- 
proach (the Dirac-Brueckner-Hartree-Fock method) and 
the other models used in our comparison (the conven- 
tional Brueckner-Hartree-Fock method together with mi- 
croscopic three-body forces) are revisited and discussed. 



II. THE "AB INITIO" APPROACH 
A. The two-body framework 

Our present knowledge of the nuclear force is the result 
of decades of struggle [J]. After the development of QCD 
and the understanding of its symmetries, chiral effective 
theories [j| became popular as a way to respect the sym- 
metries of QCD while keeping the degrees of freedom (nu- 
cleons and pions) typical of low-energy nuclear physics. 
However, chiral perturbation theory (ChPT) has definite 
limitations as far as the range of allowed momenta is con- 
cerned. For the purpose of applications in dense matter, 
where higher and higher momenta become involved with 
increasing Fermi momentum, ChPT is inappropriate. A 
relativistic, meson-theoretic model is the better choice. 

The one-boson-exchange (OBE) model has proven very 
successful in describing nucleon-nucleon (NN) data in free 
space and has a good theoretical foundation. Among 
the many available OBE potentials (some being part 
of the "high-precision generation" [1, 0, 0), we seek a 
momentum-space potential developed within a relativis- 
tic scattering equation, such as the one obtained through 
the Thompson three-dimensional reduction of the Bethe- 
Salpeter equation. Furthermore, we require a potential 
that uses the pseudovector coupling for the interaction of 
nucleons with pseudoscalar mesons. With this in mind, 
as well as the requirement of a good description of NN 
data, Bonn B Q has been our standard choice. As is 
well known, the NN potential model dependence of nu- 
clear matter predictions is not negligible. The saturation 
points obtained with different NN potentials move along 
the famous "Coester band" depending on the strength 
of the tensor force, with the weakest tensor force corre- 
sponding to the largest attraction. For the same reason 
(that is, the role of the tensor force in nuclear matter), 



FIG. 1: Contribution to the NN interaction from virtual pair 
excitation. Upward- and downward- pointing arrows represent 
nucleons and antinucleons, respectively. Dashed lines denote 
mesons. 

the potential model dependence is strongly reduced in 
pure (or nearly pure) neutron matter, due to the absence 
of isospin-zero partial waves. 

Already when QCD (and its symmetries) were un- 
known, it was observed that the contribution from the 
nucleon-antinucleon pair diagram, Fig. 1, is unreason- 
ably large when the pseudoscalar (ps) coupling is used, 
leading to very large pion-nucleon scattering lengths @. 
We recall that the Lagrangian density for pseudoscalar 
coupling of the nuclcon field (ip) with the pseudoscalar 
meson field (</>) is 

£ P s = -ig P si>^5^(f'- (1) 

On the other hand, the same contribution (Fig. 1) is 
heavily suppressed by the pseudovector (pv) coupling (a 
mechanism which became known as "pair suppression" ) . 
The reason for the suppression is the presence of the co- 
variant derivative at the pseudovector vertex, 

C pv = ^^757 M W, (2) 
TTip S 

which suppresses the vertex for low momenta and, thus, 
explains the small value of the pion-nucleon scattering 
length at threshold [9]. Considerations based on chiral 
symmetry Q can further motivate the choice of the pseu- 
dovector coupling. We will come back to this point in the 
next subsection. 

The most important aspect of the " ab initio" approach 
is that the only free parameters of the model (namely, 
the parameters of the NN potential) are determined by 
the fit to the free-space data and never readjusted in 
the medium. In other words, the model parameters are 
tightly constrained and the calculation in the medium is 
parameter free. The presence of free parameters in the 
medium would generate effects and sensitivities which 
can be very large and hard to control. 

B. The many-body framework: Brueckner theory, 
three-body forces, and relativity 

Excellent reviews of Brueckner theory have been writ- 
ten which we can refer the reader to (see [1] and references 
therein). Here, we begin by defining the contributions 
that are retained in our calculation. Those are the low- 
est order contribution to the Brueckner series (two-hole 



FIG. 2: Three-body force due to virtual pair excitation. Con- 
ventions as in the previous figure. 

lines) and the corresponding exchange diagram. With 
the G-matrix as the effective interaction, this amounts 
to including particle-particle (that is, short-range) corre- 
lations, which are absolutely essential to even approach a 
realistic description of nuclear matter properties. Three- 
nucleon correlations have been shown to be small if the 
continuous choice is adopted for the single-particle po- 
tential [13]. 

The issue of three-body forces (TBF), of course, re- 
mains to be discussed. In Fig. 2 we show a TBF orig- 
inating from virtual excitation of a nucleon-antinucleon 
pair, known as "Z-diagram" . Notice that the observa- 
tions from the previous section ensure that the corre- 
sponding diagram at the two-body level, Fig. 1, is small 
with pv coupling. At this point, it is useful to recall 
the main feature of the Dirac-Brueckner-Hartree-Fock 
(DBHF) method, as that turns out to be closely related 
to the TBF depicted in Fig. 2. In the DBHF approach, 
one describes the positive energy solutions of the Dirac 
equation in the medium as 

fEl+m*\ 1/2 ( 1 \ 

where the effective mass is given by m* = m + Us, with 
Us an attractive scalar potential. It turns out that both 
the description of a single-nucleon via Eq. (3) and the 
evaluation of the Z-diagram, Fig. 2, generate a repulsive 
effect on the energy / particle in symmetric nuclear matter 
which depends on the density approximately as 

AE ex , (4) 

and provides the saturating mechanism missing from con- 
ventional Brueckner calculations. Brown showed that the 
bulk of this effect can be obtained as a lowest order (in 
p 2 /m) relativistic correction to the single-particle prop- 
agation [Til ]. 

The approximate equivalence of the effective-mass 
spinor description and the contribution from the Z- 
diagram has a simple intuitive explanation in the obser- 
vation that Eq. (3), like any other solution of the Dirac 
equation, can be written as a combination of positive and 
negative energy solutions. On the other hand, the "nu- 
cleon" in the middle of the Z-diagram, Fig. 2, is precisely 
a superposition of positive and negative energy states. 
In summary, the DBHF method effectively takes into ac- 
count a particular class of TBF, which are crucial for 
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III. NEUTRON STARS: RESULTS AND 
DISCUSSION 



FIG. 3: Left: three-body force arising from A-isobar excita- 
tion (thick line). Right: two-meson exchange contribution to 
the NN interaction involving A-isobar excitation. 



nuclear matter saturation, but does not include other 
TBF. 

Concerning other, more popular, three-body forces, 
Fig. 3 shows the TBF that is included in essentially 
all TBF models, regardless other components; it is the 
Fujita-Miyazawa TBF [l2j- With the addition of contri- 
butions from ttN S-waves, one ends up with the well- 
known Tucson-Melbourne TBF [13l |. The microscopic 
TBF of Ref. [l5| include contributions from excitations 
of the Roper resonance (Pn isobar) as well. 

Now, if diagrams such as the one shown on the left- 
hand side of Fig. 3 are included, consistency requires 
that medium modifications at the corresponding two- 
body level are also included, that is, the diagram on the 
right-hand side of Fig. 3 should be present and properly 
medium modified. Large cancellations then take place, a 
fact that was brought up a long time ago [3| . When the 
two-body sector is handled via OBE diagrams, the two- 
pion exchange is effectively incorporated through the a 
"meson" , which cannot generate the (large) medium ef- 
fects (dispersion and Pauli blocking on A intermediate 
states) required by the consistency arguments above. 

To contrast this point of view, we will take as the other 
elements of our comparison the EoS's from the micro- 
scopic approach of Ref. [l6| . There (and in previous work 
by the authors [IB]), the Brueckner-Hartree-Fock (BHF) 
formalism is employed along with microscopic three-body 
forces. However, in Ref. [l5| the meson-exchange TBF 
are constructed applying the same parameters as used 
in the corresponding nucleon-nucleon (NN) potentials, 
which are: Argonne V18 (V18, [13]), Bonn B (BOB, 
[3), Nijmegen 93 (N93, 0). The popular (but phe- 
nomenological) Urbana TBF (UIX, [18|]) is also utilized 
in Ref. [l6| along with the V18 potential. The parameters 
which arc not specified by the NN potential are chosen 
according to independent investigations [jj| . Convenient 
parametrizations in terms of simple analytic functions are 
given in all cases for the resulting EoS's. We will refer to 
this approach, generally, as "BHF + TBF". 

In a previous work [l9[ we compared the predictions 
of the neutron skin in 208 Pb by these models and corre- 
lated them with differences in the slope of the symmetry 
energy. Model differences become larger at high-density 
and will naturally impact neutron star predictions. 



Stellar matter contains neutrons in (3 equilibrium with 
protons, electrons, and muons. Our DBHF EoS for (3- 
equilibrated matter is given in Ref. [20| . We have ap- 
plied /3-stability in the same way to the various models 
of Ref. [li| starting from the given parametrized ver- 
sions of the respective symmetric matter and neutron 
matter EoS's. At subnuclear densities, all the EoS's 
considered here are joined with the crustal equations 
of state from the work of Harrison and Wheeler [2l[ 
(for energy densities between 10 and 10 11 g cm -3 ) and 
the work of Negele and Vautherin [22j (for energy den- 
sities less than 1.7xl0 13 g cm -3 ). The composition of 
the crust is crystalline, with light |2Lj or heavy [22j 
metals and electron gas. The DBHF equation of state 
as applied in this work, including the crust, is given 
in the tables. All neutron star properties are calcu- 
lated using public software downloaded from the website 
http://www.gravity.phys.uwm. edu/rns . 

In Fig. 4, we show the mass-radius relation for a se- 
quence of static neutron stars as predicted by the various 
models listed above. All models besides DBHF share the 
same many-body approach but differ in the two-body po- 
tential and TBF employed. The resulting differences can 
be much larger than those originating from the use of 
different many-body approaches. This can be seen by 
comparing the DBHF and BOB curves, both employing 
the Bonn B interaction. Overall, the maximum masses 
range from 1.8M S (UIX) to 2.5M© (BOB). Radii are less 
sensitive to the EoS and range between 10 and 12 km for 
all models under consideration, DBHF or BHF+TBF. 
Concerning available constraints, an initial observation 
of a neutron-star-white dwarf binary system suggested a 
neutron star mass (PSR J0751+1807) of 2.1±O.2M [H. 
Such observation would imply a considerable constraint 
on the high-density behavior of the EoS. On the other 
hand, a dramatically reduced value of 1.3±0.2M Q was 
recently reported [241 ] . which does not invalidate any of 
the theoretical models under consideration. 

The model dependence is shown in Fig. 5 for the case 
of rapidly rotating stars. The 716 Hz frequency corre- 
sponds to the most rapidly rotating pulsar, PSR J1748- 
2446 (2f| , although recently an X-ray burst oscillation at 
a frequency of 1122 Hz has been reported [26] which may 
be due to the spin rate of a neutron star. As expected, 
the maximum mass and the (equatorial) radius become 
larger with increasing rotational frequency. 

Another bulk property of neutron stars is the moment 
of inertia, /. In Fig. 6, we show the moment of iner- 
tia at different rotational speeds (again, for all models), 
whereas in Fig. 7 we display the moment of inertia cor- 
responding to the maximum mass at different rotational 
frequencies. These values are not in contraddiction with 
observations of the Crab nebula luminosity. From that, a 
lower bound on the moment of inertia was inferred to be 
I >4-8 x 10 44 g cm 2 , see Ref. [27J and references therein. 
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FIG. 4: (color online) Static neutron star mass-radius relation 
for the models considered in the text. 



The size of lM mcL!C changes from model to model in line 
with the size of the maximum mass, see Fig. 7. (We re- 
call that the variations among radii from different EoS's 
are relatively mild). 

Lastly, we calculate the gravitational redshift predicted 
by each model. The redshift is defined as 



1 



2M\ 



-1/2 



1 



(5) 



This simple formula can be derived considering a pho- 
ton emitted at the surface of a neutron star and moving 
towards a detector located at large distance. The pho- 
ton frequency at the emitter (receiver) is the inverse of 
the proper time between two wave crests in the frame of 
the emitter (receiver). Assuming a static gravitational 
field, and writing goo as the metric tensor component 
at the surface of a nonrotating star yield the equation 
above. Naturally the rotation of the star modifies the 
metric, and in that case different considerations need to 
be applied which result in a frequency dependence of the 
redshift. We will not be calculating the general case here. 

From Fig. 8, it appears that the gravitational redshift 
is not very EoS-dependent (compare, for instance, the 
values at the maximum mass for each model), an indi- 
cation that the gravitational profile at the surface of the 
star is similar in all models. Thus measurements of z 
may not be the best way to discriminate among EoS's. 



IV. SUMMARY AND CONCLUSIONS 

We stress the importance of ab initio calculations to 
complement the wealth of present/future experiments 
and observations aimed at constraining the less known as- 
pects of the nuclear equation of state. In that spirit, we 
have undertaken a comparison of EoS-sensitive observ- 
ables in microscopic models, which do or do not include 
explicit TBF. In a previous work [l9j we related different 
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FIG. 5: (color online) Mass-radius relation for the models 
considered in the text and for different rotational frequencies. 



predictions of the neutron skin in 208 Pb to differences in 
the slope of the symmetry energy. An accurate measure- 
ment of the neutron skin, as expected to be taken at the 
Jefferson Laboratory (28|, is the most promising way to 
discriminate among EoS's in the low-density regime. 

At the high-densities probed by neutron stars, the 
model dependence is even larger, however, presently 
available constraints are still insufficient to discriminate 
among these EoS's. 

At high density (p between five and ten times normal 
density), the most repulsive symmetric matter energies 
are produced with BOB, V18, DBHF, N93 and UIX, in 
that order. That is consistent with the maximum mass 
predictions, which depend mostly on the absolute repul- 
sion present in the symmetric matter EoS. In pure neu- 
tron matter, N93 follows right after BOB (again, from 
largest to smallest repulsion). This indicates a some- 
what different balance when only T=l contributions are 
included. Finally, the symmetry energy, which depends 
entirely on the repulsion of neutron matter relative to 
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symmetric matter and whose density dependence con- 
trols observables such as the neutron skin, is largest in 
N93, followed by BOB, V18, UIX, and DBHF. 

The model dependence we observe comes from two 
sources, the two-body potential and the many-body ap- 
proach, specifically the presence of explicit TBF or Dirac 
effects. The dependence on the two-body potential is 
very large. Typically, the main source of model depen- 
dence among NN potentials is found in the strength of the 
tensor force. Of course, differences at the two-body level 
impact the TBF as well, whether they are microscopic or 
phenomenological (as in the case of UIX). 

On the other hand, when comparing DBHF and BOB, 
we are looking at differences stemming from the many- 
body scheme, as the two models share the same NN po- 
tential. In the BOB model, repulsion grows at a much 
faster rate than in DBHF, and more strongly so in neu- 
tron matter. (Hence, the much larger symmetry energy 
with BOB 19]). As an example, at about 6 times normal 
density the DBHF energy of symmetric matter is 67% of 
the BOB energy and only 51% at ten times normal den- 
sities. In pure neutron matter, those ratios become 49% 
and 30%, respectively. Thus, the inclusion of the micro- 
scopic TBF in the BHF model introduces considerable 
more repulsion than the Dirac effects (whose density de- 
pendence is approximately as given in Eq. (4)), through 
highly non- linear terms [161 ]. Both attractive and repul- 
sive TBF are required for a realistic description of the 
saturation point. The density dependence of the repul- 
sive terms is obviously stronger and thus dominate at 
high density. Furthermore, it appears that this is espe- 
cially true in neutron matter. The much larger (com- 
pared with DBHF) repulsion in neutron matter (relative 
to symmetric matter) seen in all BHF+TBF models may 
be due to the attractive nature of some T=0 NN partial 
waves, most noticeably Si, which, on the other hand, 
would become repulsive at high density in the presence 
of Dirac effects. This attractive component (not present 



TABLE I: Equation of state for /3-equilibrated matter, 
subnuclear densities the crustal EoS from Refs. [2ll . 
given. At normal to super-high densities the DBHF equatio 
of state is shown. 



in neutron matter) can moderate the repulsion arising 
from TBF. 

At this time, available constraints cannot pin down the 
high-density behavior of the EoS. Nevertheless, our point 
is that microscopic models allow for a deeper insight into 
the origin of the observed physical effects, and should be 
pursued along with improved constraints. 
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24200880E + 10 





12701340E4-28 





22959600E+34 





38534730E + 10 





22901680E4-28 





28966900E+34 





48625310E+10 





30733120E4-28 





36545800E+34 





61358160E + 10 





41227540E4-28 





46107400E+34 





77425110E + 10 





55287600E+28 





58169000E+34 





97699500E + 10 





74121190E+28 



TABLE I, cont. 



TABLE I, cont. 



Baryon density ( 1 / cm 3 ) 


Energy density ( g/ cm 3 ) 


Pre 


sure(d-yne/ cm^ ) 


o 


73385500E + 34 


o 


12328280E + 11 





99344590E+28 


o 


14734400E + 35 


o 


24770450E + 11 





23888140E+29 


o 


18587600E-f-35 


o 


31256660E + 11 





31991950E+29 


o 


23448300E + 35 




39441580E + 11 





42838370E+29 


o 


29580200E + 35 


o 


49769610E + 11 





57354730E+29 


o 


37313900E-+-35 


o 


62801930E + 11 





76780650E+29 


o 


47069400E-|-35 


o 


79247170E + 11 





10277560E+30 


o 


59374600E + 35 


o 


99998600E + 11 





13755860E+30 


o 


10899500E + 36 


o 


18374080E + 12 





15080040E+30 


o 


20791800E-f 36 


o 


35086530E + 12 





41327350E+30 


Q 


30685800E4-36 


o 


51813430E4-12 





66738040E+30 


o 


40581200E + 36 


o 


68549240E + 12 





89490400E+30 


o 


50477400E + 36 


o 


85290930E + 12 





10966500E+31 


o 


60373600E + 36 


o 


10203670E + 13 





12773240E+31 


o 


80166500E-|-36 


o 


13553720E + 13 





15941520E+31 


o 


90062000E + 36 


o 


15229090E + 13 





17377690E+31 


o 


10985200E + 37 


o 


18580330E + 13 





20087290E+31 


o 


11974500E4-37 


o 


20256030E + 13 





21397550E+31 


o 


13953200E + 37 


o 


23608 150E + 13 





23994040E+31 


o 


15931600E + 37 


o 


26960450E + 13 





26615840E+31 


o 


16920800E + 37 


o 


28636860E + 13 





27951740E+31 


o 


17910000E + 37 


o 


30313280E + 13 





29309910E+31 


o 


19888300E + 37 


o 


33666280E + 13 





32106350E+31 


o 


20877300E + 37 


o 


35342880E4-13 





33549110E+31 


o 


21866300E-1-37 


Q 


37019650E-+-13 





35024070E+31 


o 


23844100E + 37 


o 


40373370E + 13 





38075420E+31 


o 


24833000E + 37 


o 


42050320E + 13 





39653560E+31 


o 


25821800E + 37 


o 


43727450E-f-13 





41267750E+31 


o 


26810600E-}-37 


o 


45404400E4-13 





42918160E+31 


o 


27799300E + 37 


o 


47081710E + 13 





44605250E+31 


o 


28788000E + 37 


o 


48758840E + 13 





46329190E+31 


o 


29776600E + 37 


o 


50436140E + 13 





48089980E+31 


o 


30765200E + 37 


o 


52113630E4-13 





49887950E+31 


o 


31753700E + 37 


o 


53790930E + 13 





51722760E+31 


o 


32742200E + 37 


o 


55468420E + 13 





53594420E+31 


o 


33730600E + 37 


o 


57146080E + 13 





55502940E+31 


o 


34718900E4-37 


o 


58823750E4-13 





57447980E+31 


o 


35707300E + 37 


o 


60501410E + 13 





59429390E+31 


o 


36695600E + 37 


o 


62179070E + 13 





61447010E+31 


o 


37683800E + 37 


o 


63856910E + 13 





63500680E+31 


o 


38672000E + 37 


o 


65534760E4-13 





65589600E+31 


o 


39660200E-f-37 


o 


67212770E4-13 





67714090E+31 


o 


40648400E + 37 


o 


68890790E + 13 





69873660E+31 


o 


41636500E + 37 


o 


70568810E + 13 





72068480E+31 


o 


42624500E-f-37 


o 


72246830E + 13 





74297430E+31 


o 


43612600E-+-37 


o 


73925030E4-13 





76560670E+31 


o 


446 005 00E + 37 


o 


75603230E + 13 





78857870E+31 


o 


46576400E + 37 


o 


78959800E + 13 





83552890E+31 


o 


47564200E-f-37 


Q 


80638360E4-13 





85949750E+31 


o 


48552100E-f-37 


o 


82316 740E + 13 





88379610E+31 


o 


49539800E + 37 


o 


83995290E + 13 





90841830E+31 


o 


50527400E + 37 


o 


85673840E + 13 





93336260E+31 


o 


51515000E + 37 


o 


87352400E + 13 





95862580E+31 


o 


52502600E + 37 


o 


89031130E4-13 





98420450E+31 


o 


53490100E + 37 


o 


90709860E + 13 





10100920E+32 




54477600E-|-37 




92388600E + 13 





10362880E+32 





55465100E + 37 





94067330E + 13 





10627910E+32 





56452500E + 37 





95746240E + 13 





10896010E+32 





57439900E + 37 





97425150E + 13 





11167080E+32 





58427200E+37 





99104240E + 13 





11441110E+32 





59414500E+37 





10078320E + 14 





11718180E+32 





60401800E + 37 





10246220E + 14 





11998160E+32 





61389000E + 37 





10414150E + 14 





12280980E+32 





63363400E + 37 





10749990E + 14 





12855310E+32 





65337700E + 37 





11085860E + 14 





13440860E+32 



Baryon density ( 1 / cm 3 ) 


Energy density ( g/ cm 3 ) 


Pre 


sure(d-yne / cmP' ) 


o 


66324800E+37 


o 


11253780E + 14 





13737910E4-32 


o 


67311800E+37 


o 


11421730E + 14 





14037590E+32 


o 


69285800E+37 


o 


11757640E + 14 





14645200E4-32 


o 


70272800E+37 


o 


11925600E + 14 





14953050E+32 


o 


72246600E+37 


o 


12261520E + 14 





15576620E4-32 


o 


73233400E+37 


o 


12429500E + 14 





15892440E+32 


o 


75207000E+37 


o 


12765460E + 14 





16531580E4-32 


o 


76193800E+37 


o 


12933460E + 14 





16855220E+32 


o 


77180500E+37 


o 


13101440E + 14 





17181270E4-32 


o 


79153800E+37 


o 


13437440E + 14 





17840720E+32 


o 


80140400E+37 


o 


13605450E + 14 





18174290E4-32 


o 


821 13400E+37 


o 


13941490E + 14 





18848650E+32 


o 


83099900E+37 


o 


14109500E + 14 





19189430E+32 


o 


85072800E+37 


o 


14445550E + 14 





19878210E+32 


o 


86059200E+37 


o 


14613590E + 14 





20226200E+32 


o 


87045500E+37 


o 


14781620E + 14 





20576440E+32 


o 


88031800E+37 


o 


14949670E + 14 





20929080E+32 


o 


89018100E+37 


o 


15117700E + 14 





21284120E+32 


o 


90990300E+37 


o 


15453810E + 14 





22001090E+32 


o 


91976300E+37 


o 


15621880E + 14 





22362710E+32 


o 


93948300E+37 


o 


15958020E + 14 





23093300E+32 


o 


94934200E+37 


o 


16126090E + 14 





23461640E+32 


o 


98877600E+37 


o 


16798420E + 14 





24958070E4-32 


o 


14590200E+38 


o 


24523600E + 14 





41632200E+32 


o 


23168800E+38 


o 


38982700E + 14 





82280400E4-32 


o 


34584300E+38 


o 


582501 00E + 14 





14210400E+33 


o 


49242100E+38 


o 


83026900E + 14 





23006100E4-33 


o 


67547500E+38 


o 


11401900E + 15 





40981400E+33 


o 


78194600E+38 


o 


13207300E + 15 





54012600E4-33 


o 


89905 700E+38 


o 


15195500E + 15 





71172400E+33 


o 


10273 100E+39 


o 


17376500E + 15 





97049800E4-33 


o 


1 1672200E+39 


o 


19760500E + 15 





13704600E+34 


o 


13192900E+39 


o 


22358900E4-15 





20162700E+34 


o 


14840200E+ 39 


o 


25183700E + 15 





26192800E+34 


o 


16619200E+39 


o 


28248100E + 15 





33029700E+34 


o 


18535000E+39 


o 


31568300E + 15 





51720500E+34 


o 


20592 700E+39 


o 


35161900E + 15 





75702800E+34 


o 


22797300E+39 


o 


39049500E + 15 





11003700E+35 


o 


25153800E+39 


o 


43255500E + 15 





16084300E+35 


o 


27667400E+39 


o 


47813800E + 15 





23509200E+35 


o 


30343200E+39 


o 


52760100E + 15 





33455700E+35 


o 


33186100E+39 


o 


58135500E + 15 





46583400E4-35 


o 


36201200E+39 


o 


63986100E + 15 





63214300E+35 


o 


39393 700E+39 


o 


70359400E + 15 





83222800E4-35 


o 


42768500E+39 


o 


77303000E + 15 





10713900E+36 


o 


46330800E+39 


o 


84870300E4-15 





13477500E+36 


o 


50085600E+39 


o 


93108800E + 15 





16560600E+36 


o 


54038000E+39 


o 


10206900E + 16 





20027100E+36 


o 


58193000E+39 


o 


11181400E + 16 





24141700E+36 


o 


62555 700E+39 


o 


12243400E + 16 





29312500E+36 


o 


67131200E+39 


o 


13404700E + 16 





35869100E+36 


o 


71924500E+39 


o 


14679700E + 16 





44217300E+36 


o 


76940800E+39 


o 


16084100E + 16 





53909800E+36 


o 


82185000E+39 


o 


17629800E + 16 





64928300E+36 


o 


87662200E+39 


o 


19333700E + 16 





77966700E+36 


o 


93377600E+39 


o 


21214100E + 16 





93298700E+36 


o 


99336100E+39 


o 


23293000E + 16 





11139500E+37 





10554300E+40 


o 


25594100E+16 





13245200E+37 




11200300E+40 




28143700E + 16 





15691600E+37 





11872100E+40 





30972100E + 16 





18569700E+37 





12570300E+40 





34114400E + 16 





21904600E+37 





13295400E+40 





37606000E + 16 





25722700E+37 





14047800E+40 





41488100E + 16 





30137700E+37 





14828000E+40 





45806700E + 16 





35220200E+37 





15636600E+40 





50613000E + 16 





41065500E+37 





16474100E+40 





55963600E + 16 





47754300E+37 





17341000E+40 





61920200E + 16 





55414900E+37 





18237800E+40 





68552200E + 16 





64123300E+37 
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